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ABSTRACT 



Pulsating stars in eclipsing binary systems are powerful tools to test stellar models. Binarity enables to constrain the pulsating com- 
ponent physical parameters, whose knowledge drastically improves the input physics for asteroseismic studies. The study of stellar 
oscillations allows us, in its turn, to improve our understanding of stellar interiors and evolution. The space mission CoRoT discov- 
ered several promising objects suitable for these studies, which have been photometrically observed with unprecedented accuracy, but 
needed spectroscopic follow-up. A promising target was the relatively bright eclipsing system CoRoT 102918586, which turned out 
' to be a double-lined spectroscopic binary and showed, as well, clear evidence of y Dor type pulsations. 

• ' With the aim of combining the information from binarity and pulsation and fully exploit the potential of CoRoT photometry we ob- 

i-^ , tained phase resolved high-resolution spectroscopy with the Sandiford spectrograph at the McDonald 2.1m telescope and the FEROS 

spectrograph at the ESO 2.2m telescope. Spectroscopy yielded both the radial velocity curves and, after spectra disentangling, the 
component effective temperatures, metallicity and line-of-sight projected rotational velocities. The CoRoT light curve was analyzed 
' with an iterative procedure, devised to disentangle eclipses from pulsations. The eclipsing binary light curve analysis, combined to 

the spectroscopic results, provided an accurate determination of the system parameters, and the comparison with evolutionary models 
^ . strict constraints on the system age. Finally, the residuals obtained after subtraction of the best fitting eclipsing binary model were 

analyzed to determine the pulsator properties. 

We achieved a quite complete and consistent description of the system. The primary star pulsates with typical y Dor frequencies and 
shows a splitting in period which is consistent with high order g-mode pulsations in a star of the corresponding physical parameters. 
^ ■ The value of the splitting, in particular, is consistent with pulsations in f = 1 modes. 

■ 

. Key words. Binaries: eclipsing - Binaries: spectroscopic - Stars: oscillations - Stars: individual: CoRoT 102918586 

1. Introduction pulsations, including the non-radial pulsations typical of ^ Set, 7 

Dor and Slowly Pulsating B stars (SPB). 

^ In the last six years the CoRoT space mission has acquired more p^^^j^^^ .^^. ^^^ ^^^U^ structures can be obtained from 

,—1 than 150,000 high-accuracy light curves with dense and almost . • 1 j 1 u e ^ ,1 

' ' . ° . .-^ ° , ^ ^ , asteroseismology, and pulsating EBs have a fundamental asset: 

> ?T'!"''u''' s^"^Pli"g o^*^; t"^^ i"!^''^^!^ f™"^ fi^^ studying non-radial oscillations in EB components has the ad- 

^ While the main goals of CoRoT are exoplanet hunting by the ^^^^ ^^^^^^ ^^^^ independently derived, 

transit technique and asteroseismology of bright stars, a precious ,. -1 »u j u u- ..u ■ e 

, , . , , • , , r , with a pure geometrical method, by combining the information 

^ by-product IS a large number of vanable stai-s, most of them new ^^^^ ^j^^ j.^j^^ ^^^.^^ ^^^^^.^^ cur ves, and with un- 

- - discoveries. certai nties, in the best cases, below 1% (e.g.. ISouthworth et al.1 
EcUpsing binaries (EBs) outnumber all other variables found |2005h . Moreover, additional useful constraints derive from the 
by CoRoT and some among them reveal additional regular vari- requirement of same chemical composition and age. As the pre- 
ability superimposed on the binary light curve. The origin of the cise measurement of the mass and radius poses powerful con- 
additional variability is very often rotation modulated surface straints on the pulsational properties, EBs with pulsating compo- 
activity, in some cases - however - there is clear evidence of nent potentially provide direct tests of the modeUng of complex 

dynamical processes occurring in stellar interiors (such as mix- 

ing length, convective overshooting, diffusion). Furthermore, a 

* Based on photometry collected by the CoRoT space mission and close eclipsing binaiy containing a pulsating component is the 

spectroscopy obtained with the Sandiford spectrograph attached at the jdeal laboratory to study the effect of tides on stellar pulsations. 
2.1m telescope at McDonald Observatory, USA and the FEROS spec- 
trograph attached to the ES02.2m telescope at ESO, La Silla, Chile. The ti'ade off with these advantages is in a more complex 

The CoRoT space mission was developed and is operated by the French structure and analysis of the data, as it is necessary to disen- 

space agency CNES, with participation of ESAs RSSD and Science tangle the two phenomena at the origin of the observed time 

Programs, Austria, Belgium, Brazil, Germany and Spain series, but the results for the first CoRoT targets of these type 
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Studied so far (e.g., MaceronietalJ 120091: iDamiani et al ] l2oTot 
ISok olovskv e t al.l 2O10l) are quite encouraging. 



Most EBs observed by CoRoT belong to the exoplanet 
field, whose targets are in the V-magnitude range 1 1-16.5. Each 
CoRoT run provides the light curves of several thousands exo- 
targets (up to twelve thousands for the first runs and about half 
that number after the loss of one of the two data processing units 
in March '09). The variable stars are a utomatically classified by 
the CoRoT Variability Classifier (De bosscher et al.f 2009) which 
yields the membership probability for 29 different variability 
classes. 

The detailed inspection of the light curves of EBs and of pul- 
sating stars of the first CoRoT runs allowed us to select a small 
number of variables showing both eclipses and oscillations. All 
of them were new discoveries and the only available spectro- 
scopic measurements (if any) were low-dispersion spectra from 
a ground based survey of CoRoT fields (see next section). We 
organized, therefore, a program aimed to obtain phase resolved 
echelle spectroscopy of the suitable targets in the sample. 

In this paper we present the results relative to the brightest 
object of our sample: CoRoT 102918586, which was classified 
as EE and y Dor candidate. 

The J Dor variables are A-F stars of luminosity class IV-V 
pulsating in high-order gravity-modes, with typical frequencies 
in the range 0.3-3 days. The driving mechanism of pulsations is 
the modulation of the rad iative flux by convection at the base o f 
the convective envelope (iGuzik et alJboOOt : iDupret et al.ll2005l) . 
The g-modes probe the deep stellar interior, and - for this reason 
- are of great interest for asteroseismic studies. 

Before the CoRoT and Kepler space missions sixty-six j Dor 
variables were known, according to Henrv et al. (2007), and 
about 50 % of them were found in binary stars (12 visual bina- 
ries, 10 double-Uned and 6 single-lined spectroscopic binaries, 
no EE). The first y Dor in an eclipsing binary was announced 
by Ibanoglu et alj (1200 7') but the quality and the coverage of 
their ground-based photometry was insufficient for a detailed 
asteroseismic study. Two new ca ndid ate y Dor's in EEs were 
identified by Damiani et all ( 12010 ') and'Sokolo vskv et all (1201 Ol) 
(CoRoT 102931335 and CoRoT 102980178) but their faintness 
prevented the collection of spectroscopic data and, therefore, the 
full exploitation of the EE assets. Many others similar objects 
have been later identified in the CoRoT and Kepler archives, but 
again, without spectroscopic information. 

CoRoT 102918586 is, therefore, the first y Dor in an EE 
which can be studied in great detail, thanks to the quality of the 
CoRoT photometry, to its brightness - making possible to ac- 
quire high dispersion spectroscopy - and to the fortunate occur- 
rence of being a double-lined spectroscopic binary. 

This paper is organized as follows: Section|2]and|3]describe 
the available data and their reduction. Section |4] deals with the 
light and radial velocity curve analysis. Section |5] is devoted to 
the analysis of the disentangled component spectra, providing 
the atmospheric parameters, Section|6]presents a comparison be- 
tween the physical elements derived from the analysis and stellar 
evolutionary models and, finally. Section [T] analyzes the pulsa- 
tional properties of the primary component. 

2. CoRoT photometry 

CoRoT 102918586 is one of the ~ 10000 exoplanet targets ob- 
served by CoRoT during the "Initial Run" (IRaOl), the first 
scientific run which lasted for about 60"*. The little pre-launch 
information on this target was collected in Exo-dat database 
(operated at CeSAM/LAM, Marseille, France on behalf of the 



CoRoT/Exoplanet program, iDeleuil et al.ll2009l) . The target is 
relatively bright (V = 12.45 ± 0.03). Prior to CoRoT observa- 
tions its binary nature was unknown. We der ived a sp ectra l clas - 
sification FO V by carefully re-analyzing ( Maceroni et al.ll20 lOh 
low-resolution (R ^ 1300) spectra, which were obtained with 
the AAOmega multi-object facilities of the Anglo-Australian 
Observatory, in the frame of a program aimed t o acquire a first 
spectroscopic snap-shot of the CoRoT fields dGuenther et al.l 
2012; Sebast ian et alj|20T2l) . This spectral classification is ex- 
cellent agreement with the results obtained, later on, from the 
Sandiford and FEROS observations. 

The point-spread function (PSF) on the CoRoT exoplanet 
field detectors is actually a mini-spectrum (R~3), thanks to an 
objective prism in the optical path of the CoRoT exoplanet 
channel. Three-colour photometry (the so called red, green, and 
blue colours) is obtained by splitting the PSF into three sub- 
integrating areas for target stars brighter than V < 15.2. Those 
colours do not correspond to a ny standard photometric systems, 
though dAuvergne et al.l l2009h . Yet, they are useful to photo- 
metrically rule out false positives mimicking a planetary tran- 
sit signal (Carone et al. 2012). Therefore, even if chromatic light 
curves were available for our target, our analysis is based mainly 
on the white light curve, which was obtained summing up all the 
channels. 

This original light curve contains ~ 78000 points, with a first 
section (HJD < 2454162'') sampled in the long integration mode 
(512**) and the rest in the short one (32^*). Having checked the 
absence of high frequency components, we rebinned the time 
series to the longer step and corrected a long term trend, dividing 
the light curve by a third order Legendre polynomial. Besides, 
we discarded with a sigma-clip algorithm the obvious outliers. 
The result of these preliminary operations is shown in Fig.lT] the 
total number of points is reduced to 8062. Hereafter we will refer 
to this curve as LCO. 

The first step of our iterative procedure was the pre- 
whitening of LCO, w ith the result of the harmonic fit ob- 
tained with Period04 (lLenz&Eregedl2005h . i.e. a non-linear 
least square fit of the frequency, amplitude and phase of the 
simusoidal components, as derived from the Discrete Fourier 
Transform (DFT) of the time series. Period04 was applied to 
LCO after eclipse masking. The harmonic fit was obtained with 
the twenty-four frequencies havin g an amphtude S/N> 4, accord- 
ing to the criterion proposed by Ereger et al.l (Il993l . The first 
thirteen frequencies, in order of amplitude, are listed in Table [1] 
together with corresponding amplitudes, phases and S/N ratios. 

In this type of analysis of a single light curve a problem of- 
ten arises concerning the nature of the orbital period harmonics. 
On one side, in fact, there could be pulsations excited at orbital 
overtones by tidal action. On the other features that belong to the 
orbital variation (e.g. out of eclipse variations due to proximity 
effect or eccentricity) produce orbital overtones in the frequency 
spectrum. Besides any deviation of the binary model from the 
real configuration generates orbital overtones in the residuals 
obtained by subtraction of the binary model. The analysis of 
CoRoT 102918586 is a good example in this respect: the first 
inspection of the light curve, and the knowledge about the lumi- 
nosity class of the primary star, suggested a binary configuration 
with two components of similar surface brightness (and hence 
almost equally deep minima) in a circular orbit, and with an or- 
bital period of ^ 8'?78. This implies, as well, small fractional 
radii of the components and very little tidal deformation. 

While it was possible t o obtain an excellent fi t of binary + 
pulsations with this model dMaceroni et al.ll2"010l) . it was prob- 
lematic to explain the presence of several orbital overtones in 



2 



C. Maceroni et al.: CoRoT 102918586: a y Dor pulsator in a short period eccentric eclipsing binary. 



the harmonic fit (see Table [TJ, as tidal forces should be weak in 
spherical components describing circular orbits. 
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Fig.l. CoRoT 102918586 de-ti-ended white-light Hght-curve 
after rebinning to 512^ and normalization to the mean value 
(echpses excluded). 

The alternative hypothesis, suggested in the above men- 
tioned paper, was a completely different system configuration: 
the minima seen in the curve are the same, single eclipse. In 
that case the system should be formed either by very diflFerent 
stars, or have an elliptic orbit with such an orientation in space 
to produce an eclipse only when the components are close to 
periastron passage (e .g., a case similar to that of HD 174884, 
iMaceroni et al.ll20()9h . This second hypothesis resulted to be the 
true one after collection of the first series of spectra: the orbit 
turned out to be elliptic, with half the orbital period, and its incli- 
nation angle and longitude of periastron yield a single, grazing, 
eclipse (see Section|3]l. 



Table 1. Results of the preliminary harmonic analysis. 



F(d-') 


Ampl -10' 


Phase (2;r rad) 


S/N° 


remark 


1.22452 (2) 


4.30116 (1) 


0.1576 (3) 


24.7 




1.12666 (3) 


3.05882 (1) 


0.1223 (4) 


21.7 




1.17256 (3) 


2.64945 (1) 


0.6155 (5) 


24.2 




0.94653 (4) 


1.43687 (1) 


0.2809 (9) 


16.5 




2.3512(1) 


0.72474 (1) 


0.101 (2) 


11.8 




0.4572 (1) 


0.66486 (1) 


0.432 (2) 


8.51 


2 Foi-b 


0.2277 (1) 


0.53671 (1) 


0.286 (2) 


7.00 


Forb 


2.39622 (1) 


0.49724 (1) 


0.903 (3) 


9.77 




0.0519(2) 


0.50943 (1) 


0.251 (3) 


6.72 




2.4490 (2) 


0.41644 (1) 


0.097 (3) 


9.65 




0.6833 (2) 


0.44263 (1) 


0.055 (3) 


7.96 


3 Foi-b 


2.3009 (2) 


0.42179 (1) 


0.944 (3) 


9.33 




0.9118 (2) 


0.32531 (1) 


0.901 (4) 


6.49 


4F„,b 



Number in parenthesis: formal errors on the last digit of the LS fit. 
" S/N is computed over an interval of 5 d"' 



Fig. 2. Upper panel: the phased light curve of the eclipsing bi- 
nary as residuals of the harmonic fit, orbital overtones excluded. 
Full line: final fit with PHOEBE. Lower panel: fit residuals. The 
displayed curve is the final result of our procedure (after two 
steps of iterative pre-whitenening, see text). 

the sides of the single minimum, which is exactly the expected 
shape for a longitude of periastron close to 90 degrees, i.e. a line 
of view coplanar with the system semi-axis. 

The light curve of Fig. |2]is phased with the ephemeris: 

r^ini = 2454139.3798(4) -i- 4^39138(6) x E. (1) 

which was derived by a linear fit to the minimum times, com- 
puted in their turn with a parabolic fit to the lower section 
of each minimum. The figure also shows the synthetic light 
curve of the best fittin g binary model computed with PHOEBE 
(iPrsa & Zwittedl2005h and described in Section |4] 

For the sake of brevity we directly show the final result of the 
iterative process used to extract LCEB from LCO. The procedure 
yielding the final LCEB is: 

1 . prewhitening of LCO with 20 out of the 24 frequencies ob- 
tained in the first harmonic fit, as akeady described 

2. phasing with the ephemeris of Eq.[T]and further binning with 
a variable step, to limit computing time for PHOEBE fitting 
(namely computation of normal points in phase bins of 0.001 
and 0.002 in and out of eclipse). 

3. subtraction of the PHOEBE fit model from LCO and new 
harmonic fit of this full time series. 

4. prewhitening of LCO with the results of the previous step, 
yielding a new 'orbit-only' light curve. 

The procedure was stopped after this last stage because there 
was no significant change in the parameters of the binary model 
and in the results of the harmonic fit. 

The resulting LCEB, which is phased and binned as in step 
2, is shown in Fig.|2] together with the residuals of the PHOEBE 
fit. Note that the larger scatter of residuals at primary minimum, 
and at the predicted position of the secondary one, is partly due 
to the smaller number of points in the variable width bins. 



The contribution to LCO only due to eclipses (i. e. the resid- 
uals of subtraction of the harmonic fit), was then recomputed 
subtracting only the 20 frequencies not corresponding to or- 
bital overtones. This orbit-only light curve (Fig. |2] LCEB) is no 
longer flat out of eclipse and shows two almost equal bumps on 



3. High resolution spectroscopy of 
CoRoT 102918586 

CoRoT 102918586 is the brightest target of the sample of 
CoRoT binaries with pulsating components we selected for spec- 
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troscopic follow-up. We collected echelle spectra with two dif- 
ferent instruments, namely: the Sandiford spectrograph, attached 
on the Cassegrain focus of the 2.1m Otto Struve telescope at 
McDonald Observatory (Texas, USA) and the fiber-fed FEROS 
spectrograph mounted on the ESO 2.2m telescope in La Silla, 
Chile. 

Seven Sandiford spectra were acquired in December 2009 
and three in January 2011. We used a 1 " slit, which gives a 
resolving power of = 47 000. The small size of the CCD de- 
tector enables to cover only ~ 1000 A. We thus set the spectro- 
graph's grating angle to encompass the 5050 - 6030 A wave- 
length spectral region, which is a good compromise between ob- 
tained signal-to-noise (S/N) ratio and number of spectral lines 
suitable for radial velocity (RV) measurements. Two consecu- 
tive spectra of 20-30 minutes were secured in each observing 
night and subsequently co-added to remove cosmic rays hits. 
The inspection of the first seven spectra, acquired at McDonald 
in December 2009, revealed that CoRoT 102918586 is a double- 
lined spectroscopic binary with an orbital period of about 4.39 
days. 

Seven more spectra were taken with FEROS in December 
2010. FEROS is fed with a 2 " fibre, which yields a resolution 
of 7? = 48 000 and a spectral coverage from 3800 to 8500 A. We 
adopted the same observing strategy as that used in the Sandiford 
obsei-vations. 

The reduction process involved the usual steps of de-biasing, 
flat-fielding, background subtraction and wavelength calibration 
by means of measurements of a ThAr calibration lamp. The 
spectral section of interest were then normalized to the contin- 
uum and baiycentiic corrections were computed for the time of 
mid-exposure. The FEROS spectra had a typical S/N of 50-60 in 
the region 5500-5800 A andof 30-40at4800 A, while Sandiford 
spectra have a typical S/N 40-50 at 5500 A. 

The component radial velocities were derived with a stan- 
dard cross correlation algorithm, IRAF's tool FXCOR. In the 
case of FEROS spectra, we used the wavelength interval 4500 
- 4800 A, rich of metallic lines. As cross correlation template, 
in the FEROS case, we chose a nearby bright CoRoT primary 
target (HD 49933) of similar spectra type (F2) and we used, 
for reference, the RV standard HD 22484. The cross correlation 
template used for the Sandi ford spectra wa s instead the spectro- 
scopic standard HD 50692 dNidever et al.li2002) . The radial ve- 
locities of the components were derived by fitting two Gaussian 
curves (with the FXCOR "deblend" function), obviously the fit 
eiTors depend on the separation of the CCF peaks. 

The radial velocity curves, showing excellent agreement be- 
tween the two data sets, are shown in Fig. [3] the corresponding 
values in Table |2l Both Sandiford and FEROS un certainties in- 
clude those of the RV stand ard stars (0. 1km s~^ : i Nidever et al.l 
[2002; Holmb erg etakl |200l) , the FEROS data (the last seven 
rows in Table |2]l include as well the uncertainty on the template 
velocity (0.1 kms"'). 

4. Light and radial velocity curve analysis 

4. 1. Determination of the best fitting binary model. 

The light and radial velocity curve solutions were performed 
with the current ("devel") version of PHOEBE, which includes 
flux computation with the CoRoT transmission functions, for 
both the seismology and the exoplanet field (for details see 
iMaceroni et al. 2009). 

We preferred to avoid the simultaneous solution of light and 
radial velocity curves, being the RV data acquired at different 



Table 2. Radial velocities of CoRoT 102918586 



BJD -2450000 






Vrad,2 


Instmment" 




(kms ') 


(kms-') 




5173.91504 


81.6 ±0.4 


-4.5 


±0.5 


S 


5174.92337 


128.4 ±0.3 


-51.9 


±0.3 


S 


5175.79674 


13.3 ±0.4 


79.9 ±0.4 


S 


5176.85915 


-32.9 ±0.2 


125.5 


±0.3 


S 


5177.99141 


57.2 ±1.0 


23.3 


±1.0 


S 


5178.97184 


122.5 ±0.3 


-47.9 


±0.3 


S 


5179.82042 


85.3 ±0.3 


-7.5 


±0.4 


s 


5584.90412 


-51.6 ±0.3 


143.3 


±0.3 


s 


5585.90005 


21.5 ±1.0 


63.1 


±1.0 


s 


5565.593/3 


-42.1 ±0.3 


136.1 


±0.3 


c 

a 


5556.80513 


127.3 ±0.5 


-54.2 


±1.0 


F 


5558.57947 


-49.7 ± 0.6 


144.5 


±1.2 


F 


5558.64812 


-47.5 ± 0.5 


140.3 


±1.1 


F 


5560.74155 


110.0 ±0.6 


-34.4 


±1.1 


F 


5560.76300 


111.0 ±0.5 


-36.0 


±1.2 


F 


5561.56488 


118.8 ± 0.5 


-46.1 


±1.0 


F 


5562.65659 


-50.2 ± 0.6 


141.7 


±1.1 


F 



" S for Sandiford spectrograph and F for FEROS. 
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Fig. 3. Upper panel: the phased radial velocity curves of 
CoRoT 102918586 components. Filled circles: FEROS mea- 
surements, filled squares Sandiford data; the uncertainties are 
smaller than the symbol size. Lines: PHOEBE fits (full: primary, 
broken: secondary component). Lower panel: fit residuals, with 
the same symbols. 



epochs with respect to photometry and being the data sets very 
different in terms of obsei-ved point number and accuracy. On the 
other hand the two solutions were connected by keeping fixed in 
each of them the parameters better determined by the other type 
of data. 

In the light curve solution we adjusted the inclination /, the 
secondary effective temperature, Teff^2, and the surface potentials 
Qi_2; the primary passband luminosity Li was separately com- 
puted rather than adjusted, as this allows for a smoother conver- 
gence to the minimum. The eccentricity e, mass ratio q, and the 
longitude of periastron oj were instead fixed to the values derived 
from the fit of the radial velocity curves (where we adjusted as 
well the system semi-axis, a, and the barycentric velocity y). 

For limb darkening we adopted a square root law that em- 
ploys two coefficients, xld and yio per star and per passband. 
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Table 3. Parameters from PHOEBE fits of light and RV curves 
and derived physical parameters of CoRoT 102918586 



System 

Primary Secondary 

To 77.66 ± 0.07 

e 0.249 ± 0.005 

w 102.6° ± 1.4 

q 0.898 ± 0.007 

a (Rq) 16.53 ± 0.07 

7(kms"') 42.0 ±0.4 

(i2/i'i)coRoT 0.673 ± 0.015 

Teff (K) 7400° + 90 7144± 1 10 

M(Mo) 1.66 ±0.02 1.49 ±0.03 

R{Re) 1.64 ±0.01 1.48 ±0.01 

logg 4.23 ±0.01 4.27 ± 0.01 



" Fixed value, the uncertainty is from the analysis of Section [STT] 



The coefficients are determined by PHOEBE interpolating, for 
the given atmospheric parameters, limb darkening coefficient 
tables. These tables were computed by synthesizing SEDs for 
many passbands and at different emergent angle, and fitting the 
linear cosine law and non-linear lo g and square root laws by least 
squares (see iMaceroni et al.l 12009,1 . The gravity darkening and 
albedo coefficients were kept fixed at thei r theoretical valu es, 
g - 1.0 and A = 1 .0 for radiative envelopes (Ivon Zeipell[T92l . 

The exoplanet fields were selected with the purpose of max- 
imizing the number of late type dwarfs, which, thanks to higher 
contrast, are better suited for transiting planet detection. As 
a consequence the exo-fields are usually relatively crowded. 
Moreover the target PSFs, as we have akeady mentioned, are en- 
larged by the miniprism. To limit contamination each target flux 
is recorded only in a specific (software) mask, whose size and 
form depends on the star characteristics. Nevertheless contami- 
nation by faint unresolved targets is rather common. Information 
on the degree of contamination has been independently acquired 
by ground-based photometry programs with a much higher space 
resolution (Deeg et al. 2009). For our system the contamination 
is non-negligible, being estimated to 15% of the total flux (Deeg 
2012, private communication), a value similar to that appearing 
in the Exodat database (14.4 %). We added, therefore, the former 
value as "third light". 

The primary effective temperature was kept fixed, as it is well 
known that the solution of a single light curve is sensitive only 
to relative values of the effective temperatures. We started with a 
value of 7200 K (an estimate from the spectral type) and finally 
used 7400 K (suggested later on by the analysis of the disentan- 
gled spectra). As expected the change of Tefr,i reflected into a 
similar shift for Teff.2. The analysis of the line profiles (see sec- 
tion |5]l suggested rotational velocities very close to spin orbit 
synchronization for both components, we adopted, therefore, a 
ratio of rotational to orbital period /i,2=l- 

The system model corresponding to the best fit is formed by 
two almost twin stars; the low inclination, the eccentricity and 
the orientation of the orbit are at the origin of the single mini- 
mum, which is a grazing eclipse of the larger star. The system 
parameters from the light and radial velocity curve solution are 
collected in Table |3] The unicity of the solution and the deriva- 
tion of parameter uncertainties are discussed in the next section. 



4.2. Unicity of tine soiution and parameter uncertainties 

The light curve fit is based on the minimization procedure of a 
cost function measuring the deviation between model and obser- 
vations in the space of adjusted parameters ( in PHOEBE differ- 
ential corrections in the Levenberg-Marquart's variant or Nelder 
and Mead's downhill simplex). As a consequence it suffers of the 
well know problems of these methods: the minimization algo- 
rithm can be trapped in a local minimum or degeneracy among 
parameters and data noise can transform the minimum into a 
large and flat bottomed region, or an elongated flat valley, rather 
than a single point. Besides, the correlation among the param- 
eters implies uncertainties on the derived values which are sig- 
nificantly larger than the formal errors derived, for instance, by 
the least square minimization. In the specific case of a detached 
binary with almost grazing eclipses one expects, for instance, a 
strong degeneracy among fractional radii (i.e. surface potentials) 
and inclination. 

To handle these problems we performed an heuristic scan of 
the hyper-surface describing the cost function value in the space 
of the adjusted parameters. Minimizations were first performed 
starting from random assigned input points, allowing to locate 
the region(s) of minima. The lowest minimum region was then 
further explored by mapping the cost function at randomly cho- 
sen points. Fig.|4]displays the lower envelope of the projections 
on different parameter plans of the cost-function hyper-surface. 

We chose, as uncertainty on the parameters, the interval cor- 
responding to a variation of 1% of the absolute minimum value. 
The choice of this threshold value is based on the comparison 
with the results of an independent estimate of uncertainties from 
bootstrap resampling (BR), a very useful technique to estimate 
parameter c onfidence levels o f the least squares solutions (see, 
for instance. i Press et anil992l) . In short, BR consists in generat- 
ing many different data samples by random resampling with rep- 
etitions (bootstrapping) the available data, performing the mini- 
mization procedure for each sample, and deriving confidence in- 
tervals from the resulting distribution of parameters. The interval 
of a given parameter containing 68.3% of the solutions has the 
same meaning of the l-cr interval of the Gaussian distribution. 
The estimate by BR was performed both for the radial and the 
li ght curve solution, according to th e scheme described in detail 
in IMaceroni & Rucinskil (Il997h and IMaceroni et al" ] d2009h .The 
main point is that the procedure is performed within the min- 
imum already established by a single iterated solution (that is, 
using only one set of residuals and parameter derivatives). The 
68.3% interval for the inclination from BR is +0.07. A compar- 
ison the first panel of Fig. fallows to derive the threshold value 
of 1% used in the following. 

To illustrate the correlation between the sum of radii and the 
inclination we use the quantity: 



D.i-q Q2 - 0.5 + 0.5^ 

where Qi, Q2 are the surface potentials, rs is a good approxima- 
tion of the sum of the star fractional radii in the case of spherical 
components. In Figure |4] we show the relation between rs and / 
for the sample of points with a cost function value < 1 % above 
the minimum. As expected the correlation is very strong. The 
projection on the y-axis of the distribution provides an estimate 
of the uncertainty on rs (~ 0.001) taking the correlation with 
/ into account. This is larger than the formal error on the same 
quantity, which can be derived from the LS fit results, by a factor 
of ~ 8. 
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Fig. 4. The projection of the cost function along inclination (top left), primary surface potential (top right) and secondary temperature 
(bottom, left). The minimum value is indicated by a filled dot. Bottom right panel; the correlation of rs with inclination for the best 
fitting models (cost function value < 1 % above the minimum). 



The uncertainties on the adjusted and derived parameters re- 
ported in Table [3] correspond to the 1 % threshold (and to the 
similar BR results), all significantly larger than the formal errors 
of the LS fits. The formal errors from the LS fit are typically 
a factor of ten smaller. The eiTor on the primary temperature, 
which is fixed in the solution, derives from the analysis of the 
disentangled spectra. 



5. Spectroscopic analysis 

5.1. The effective temperatures 

Due to Doppler shifts the spectra of binary stars are complex 
and vary in the course of the orbital cycle. Usually the spec- 
tral lines of the components are overlapping, especially around 
eclipses, what in turn makes the radial velocities difficult to mea- 
sure and inaccurate. The method of spec tral disentangling, (spd , 
for short), ove rcomes these difficulties (ISimon & Sturml[l992t 
lHadravalfT995h . In a self-consistent way it returns the optimal 
set of orbital elements and the individual spectra of the com- 
ponents. The gain of spd is obvious: separated spectra could 
be analyzed as single star spectra, without the interfering sig- 
nal of the companion. This allows the spectroscopic determi- 
nation of the effective temperatures and of detailed elemen- 
tal abundances (Hensberge et al. 2000; Pavlovski & Hensberge 
2005; Frematetal. 2005; Clausen et al. 2008; Hareteretal. 
2008: .Pavlovski & Southworth i2009i: iTkachenko et al.i ,2009.) . 



The effective temperature determined from hydrogen line pro- 
files, and/or from the strength of temperature sensitive lines, su- 
persede in accuracy and precision those derived from e.g. broad- 
band or Stromgren photometry. The same is true for the metal - 
licity derived from a detailed abundance study rather than from 
photometric colour indices. 

For the purpose of spectroscopic diagnostics of the compo- 
nent spectra spd was performed on the time-series of the FEROS 
spectra, since only these spectra cover the hydrogen Balmer lines 
on which is based the deter mination of the effective temper- 
atures. We used FDBinar\{3 (lliiic et an i2004l) which performs 
spectral dis entangl ing in the Fourier space according to the pre- 
scription of iHadrav a (1995). In spite of the very limited num- 
ber of spectra spd performed qui te well and, thanks to f air phase 
distribution of the obsei-vations (iHensberge et al.ll2008l) . did not 
produce strong undulations in disentangled spectra. 

The determination of the effective temperature from the pro- 
files of the hydrogen Balmer lines suffers from degeneracy in 
Teff and log g, and that is a major limitation in setting-up the 
appropriate model atmosphere for a detailed abundance study. 
However, the complementary analysis of the light and RVs 
curves of eclipsing double-lined spectroscopic binaries provides 
as well an independent and accurate determination of the surface 
gravity of the components. In well-determined solutions the ac- 



http:sail.hr.fdbinary.html 
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4850 4860 4870 4880 

wavelength 

Fig. 5. The best fitting synthetic composite spectrum of CoRoT 
102918586 (solid line in dark grey). The spectral window is cen- 
tered on the HyS line. The FEROS spectrum observed on HJD 
2455558.64812 is shown in light grey. The primary lines are 
blue-warded with respect to the secondary ones. 



curacy in log g is around 0.01 dex or less, and that is indeed our 
case, so that the degeneracy in Tetf and log g could be lifted. 

The components' spectra are diluted in a binary star spec- 
trum proportionally to their relative light contribution to the 
total light of the system. In the observed spectra of a binary 
system there is no information on the absolute spectral line 
strengths and, unless there is at least one spectrum taken in 
eclipse, there is an ambiguity in disentangled component spec- 
tra. Therefore, a renormalisation of the disentangled spectra of 
the individual components to their own continuum shall be per- 
fo rmed using some additional info rmat ion, as discussed in detail 
by Pavlo vski & Hensbergd (1201 Oh and iPavlovski & SouthworthI 
(2012). 

In the pr esent analysis we u sed a constrained optimal fit- 
ting method ( Tamaio et al.ll201 lb . In the new release of the com- 
puter code for constrained optimal fitting of disentangled com- 
ponent spectra the following parameters can be adjusted for both 
stars: effective temperature, surface gravity, light dilution factor, 
Doppler shift, projected rotational velocity, and the vertical shift 
to adjust for the continuum (in spectral disentangling the dis- 
entangled component spectra are shifted and an additi ve factor 
shal l be appli ed to return the m to the c ontinuum level (j l liiic et al.l 
12004: Pavlovski & Hensberge 2005: Tkachenko et al ll2009t) . 

The disentangled spectra were then analyzed to derive the 
atmospheric parameters of each component. To this purpose 
a comprehensive grid of theoretical spectra was calculated 
for^eff = 5000 to 15000 K, and log g = 2.5 to 5.0 in steps of 
250 K in temperature an d 0.5 i n log g . The atlas9 model atmo- 
spheres of Castelli et al .l (1 19971) and the uclsyn spectral synthesis 
package (Smallev et al. 1997) were used for the construction of 
the grid ( which is for fixed solar metallicity but was used only 
in the initial determination of Teff from t he Bakner lines, see be- 
low). A genetic algorithm was then used dCharbonneaul 1 9951) for 




CD 




5195 5200 5205 5210 

wavelength 

Fig. 6. Short spectral stretch of two observed composite spec- 
tra at different orbital phases from the FEROS dataset. The two 
individual disentangled spectra are shown at the bottom along 
the best fitting synthetic spectra calculated for the effective tem- 
peratures and surface gravities derived in Sect. 0& 2](Tab. |3]& 

a. 



the optimization with respect to model parameters. Further de- 
tailed description of the code is given in Pavlovski et al. (2012, 
in preparation). 

Only the hydrogen Balmer lines Hy and Hy6 were selected 
for the optimal fitting, excluding from the fit the blends due to 
the various metal lines. The surface gravities and the light fac- 
tors of both components were fixed according to the light curve 
solution (Table |3]l. Since the Balmer lines are almost insensi- 
tive to rotational broadening, with the exception of the very core 
of the lines, the projected rotational velocities were also kept 
fixed. Their values were derived from the clean metal lines, in 
an iterative cycle with the determination of the effective tempera- 
tures from Balmer lines yielding Vrot sin / = 20.5 + 0.5 km s"' and 
v'rot sin / - 14.5 + 0.5 km s^'for the primary and secondary com- 
ponent, respectively. With the surface gravities, the light dilution 
factors and projected rotational velocities for both components 
fixed, the only free parameters in the the optimal fitting were the 
effective temperatures, the Doppler shifts of disentangled spectra 
to reference frame of the synthetic spectra, and the vertical ad- 
justment of the continuum level for both components separately. 

The mean values of the effective temperatures derived in 
these calculations, separately for two Balmer lines read: Teff = 
7 400+90 K for the primary, and T^ff = 7 100+1 10 K for the sec- 
ondary component. The latter is in excellent agreement with the 
value derived from the light curve solution in Table |3] The qual- 
ity of the fit for the composite H/3 line in one of the observed 
spectra is shown in Fig.|5] 

5.2. Abundances 

The important outcome of spectral disentangling is the gain in 
S/N ratio. That is evident from the short section of the compo- 
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nent disentangled spectra which are displayed in the lower part 
of Fig. |6l together with two composite observed spectra at dif- 
ferent orbital phases. 

We estimated an enhancement in S/N, thanks to spectral dis- 
entangling and binning of disentangled spectra, by the factor of 
about 2.2 and 1.6, for the primary and the secondary respec- 
tively. The difference is due to the unequal contribution of the 
components to the total light. The increase in S/N makes, in turn, 
possible the abundance study. 

First, we estimated the micro turbulent velocity In the 
spectra of F-type stars Fe i lines serve to the purpose, since they 
are the most numerous. The mictroturbulent velocities for the 
primary and secondary component, minimizing the scatter in the 
abundances from different iron lines, are ^t,p = 1.8 + 0.1 kms ' 
andft,s = 1.0 + 0.1 kms"'. 

Iron appears in two ionization stages, Fe i and Fe ii, which 
were used as additio nal diagnostics of the effective tempera- 
tures (ISmallevll200 5l). The null-dependence of the iron abun- 
dance on the excitation potential slightly corrected the T^s 's, 
in a larger amount for the secondary component. However, the 
Fe I lines outnumber the Fe ii lines, and the solution is not well 
constrained. Differences of about 0.20 dex are found in the 
iron abundances from Fei and Fen, for the primary and sec- 
ondary component, which are larger than the estimated errors in 
abundance determination. In both cases the iron abundances are 
larger from the Fe ii lines. Similar differences were found also for 
other ions which appear in two ionization stages. However, the 
Fe I lines outnumber the Fe ii ones by factor of 5 and 6 for the pri- 
mary and secondary component, respectively, and we eventually 
adopted abundances derived from Fe i lines as final. Still, this en- 
ables a slight correction of Tetj 's derived in the previous subsec- 
tion, and the final, adopted Tefj 's are listed in Table |4] The final 
adopted elemental abundances, as derived in LTE approxima- 
tion, are listed in Table |5] The abundances relative to hydrogen 
are given in the second and fourth columns (log A^(H) = 12.00), 
those relative to the solar one in th e third and fifth columns , 
with the solar abundances as given bv lGrevesse & Noelsl(ll993h . 
We conclude that within Icr errors the abundances for both 
components are close to the solar one. The Iron abundances 
[Fe/H]p = 0.11+0.05 and [Fe/H]s = -0.10+0.04 give an approx- 
imate metallicity for the components, Zp = 0.020 ± 0.001 and 
= 0.016 ± 0.001, being Zq = 0.02 from iGrevesse & Noelsl 



Table 5. Photospheric elemental abundances of 
CoRoT 102918586. 



( 1199 3*). Because of the uncertainties in the light ratio between 
the components, which in turn could affect our Tetf 's determi- 
nation, we conclude that the metallicity of both components can 
be assumed to be solar 



6. Physical properties of CoRoT 102918586 

A comparison with theoretical models of the given mass and 
chemical composition provides further information on the phys- 
ical properties of the system. For instance, the inspection of the 



Table 4. Atmospheric parameters of CoRoT 102918586 
components^. 



Ion 




[X/H]p 




[X/H]s 


C 






8.51+0.09 


-0.04+0.10 


Si 


7.43±0.05 


-0.12+0.07 


7.48+0.06 


-0.07+0.08 


Ca 


6.22±0.15 


-0.14+0.15 


6.06+0.14 


-0.30+0.14 


Ti 


5.11+0.17 


0.09+0.18 


4.94+0.09 


-0.08+0.11 


Cr 


5.83+0.08 


0.16+0.09 


5.54+0.12 


-0.13+0.12 


Mn 


5.55+0.11 


0.16+0.11 


5.54+0.09 


0.15+0.09 


Fe 


7.61+0.03 


0.11+0.05 


7.40+0.03 


-0.10+0.05 


Ni 


6.37+0.08 


0.12+0.09 


6.23+0.07 


-0.02+0.08 



evolution of the stellar radii with time, combined to the con- 
straint of coevality, allows an accurate determination of the sys- 
tem age. We prefer to use the evolution of the radius, instead, 
for instance, of the position in the Hertzsprung-Russell (H-R) 
diagram, because the radii are directly and precisely determined 
from the radial velocity and light curve solutions, while the com- 
putation of temperature and luminosity implies the use of color 
transformations and bolometric corrections. 
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Teff 's are derived from the Balmer lines, v sin / and from metallic 
lines. 



Fig. 7. Evolution of the component radii. Upper panel: based on 
models without overshooting, aov = 0, lower panel: ckov = 0.2 . 
The light-gray bands describe the time evolution of the radii for 
star in the mass range Mi 2 ± AMi 2, from Table [3] The darker 
areas fill the intersection with the R\ 2±l^\,i lines from the same 
Table. The vertical lines bound the time interval respecting the 
constraint of coevality. 
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Fig. 8. The location of CoRoT 102918586 (starred symbol) in 
the period - eccentricity diagram for stars of spectral type A8-F2 
(derived from t he database of the 9*^ Catalogue of Spectroscopic 
Binary Orbits, iPourbaix et"ani2004 



The evolutionary tracks shown in Fig.|7]were obtained from 
stellar evolution modeling with the code CLES (Code Liegeois 
d'Evolution Stellaire. IScuflaire et al.ir2008h . The computations 
were made with the equa tion of state from OPAL (OPAL05, 



[Rogers & Nayfonov 2002). The op acity tables are, as well, 
from OPAL ( 'iglesias & Rogers 1996) for the solar mixture of 
iGrevesse & Noels ( 1993 , GN93), extended to low temperatures 
with the Ferguso n et al.l (l2005b opacity values. The nuclear re- 
action rates are those of the NACR E compilatio n ('Angul o|l999h 
except for '''N(p, r)'^0, updated bv lFormicola et al. (200^^ 

Models with and without overshooting were comput ed with 
the m ixing-length theory (MLT) of convection CBohm-Vitensd 
Il958l) . The value of q-mlt was kept fixed, adopting the solar 
value of 1 .8 (at any rate, in the temperature range of interest here, 
its value does not affect the stellar radius). For the models with 
overshooting, the overshoot length was expressed in terms of the 
local pressure scale height Hp, as aovHp, and the aov value was 
fixed to aov = 0.2. The chemical composition was assumed to 
be solar, according to the results of the previous Section. The 
computations were done without including microscopic diffu- 
sion. 

Figure Q shows, for both components, the region of the 
R (log t) diagram bounded by the evolutionary tracks of mass 
M ± AM from Table [3] The intersection with the lines R + AR 
of the same Table provides the possible age range for each star, 
the further constraint of coevality narrows the range to the final 
value of 540 + 120 Myr, in the case of no-overshooting, and to 
580 +100 Myr for cqv = 0.2. In both cases a young system with 
two Main Sequence (MS) stars which have burnt less than 20% 
of their initial hydrogen content. 

Because of the young age the binary had the time reach only 
spin-orbit synchronization while orbit circularization is still on 
the way. The eccentricity is quite high for the orbital period, as 
can be seen (Fig.[8]l, from the location of CoRoT 102918586 in 
the period-eccentricity diagram of stars of similar spectral type. 
For periods smaller than -10 days the circularization process is 
expected to be efficient enough to change the system eccentric- 
ity in a few Gyr. So in Fig. |8]the younger stars born with high 
eccentricity trace the upper envelope of the distribution, an will 
evolve with time towards lower and lower values. 



7. Pulsational properties 

The final frequency spectrum, after subtraction of the binary 
model, is shown in Fig. |9] The plot is restricted to the frequency 
interval containing meaningful features as, for higher frequen- 
cies, we only detected signal related to the satellite orbit. 

As mentioned before, to evaluate the frequency significance 
we used the criterion suggested by Breger et al. (1993), i.e. aS/N 
value of the amplitude > 4. The twenty-eight frequencies satis- 
fying this constraint are listed in Table [6 ] The e rrors in the table 
are derived according to lKallinger et alT(l2008l) and are sensibly 
larger that the formal errors of the LS fit, the remark column lists 
the closest frequency combinations, taking the uncertainties into 
account. 

The first three dominant frequencies of the amplitude spec- 
trum are separated by ~ 0.05 d The fourth frequency is the 
beating of the third highest amplitude frequency (F3 of Table |6]l 
and Forb = 0.22772 d Other combinations of these four 
highest amplitude frequencies are also clearly visible. Figure |9] 
shows as well the window DFT, whose main features are the 
peaks at 13.972 d^'and its day aliases, due to the satellite or- 
bital frequency, and a peak at 2.005 d 'due to the South Atlantic 
Anomaly crossing. The gray line is the significance threshold 
which is four times the local mean of the residual spectrum. 

The Foib overtones present in the initial analysis of Table [T] 
have disappeared or very weak, which is an indirect confirmation 
of the successful removal of the orbit-only binary light curve. 
The low frequency components (F < 0.1 d are presumably 
leftovers of the detrending procedure with a low order polyno- 
mial. 



0.005 




Fig. 9. The final amplitude spectrum of CoRoT 102918586 and, 
in the inset, the window DFT. The vertical dotted lines corre- 
spond to the orbital period overtones, the broken lines to 2F3 
and 3F3. The gray line is the threshold of significance (S/N> 4) 
computed as explained in the text. 



An oscillatory pattern of the kind superimposed on our EB 
light curve is not necessarily due to pulsations. A study of the 
large sample of va riables detected by the Kepler satellite and 
classified as y Dor (iBalona et al.ll201 ll) has shown that the vari- 
ability of "symmetric" type, i.e symmetric with respect to its 
mean value, may be due to stellar spots experiencing differential 
rotation, and that the Kepler sample of y Dor is contaminated by 
non-pulsating star. In this case the dominant frequency is that of 
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Table 6. The frequencies derived from the analysis of the final set of residuals, after subtraction of the EB model, with an amplitude 
S/N> 4. 





F(d-') 


Ampl -10' 


Phase (27r) 


S/N 


remark " 


Fl 


1.2247 + 0.0017 


4.25 + 0.41 


0.59 


+ 0.04 


17.8 


g 


F2 


1.1255 + 0.0012 


3.10 + 0.22 


0.16 


+ 0.03 


24.5 


g 


F3 


1.1716 + 0.0010 


2.73 + 0.16 


0.03 


+ 0.03 


29.3 


g 


F4 


0.9462 ± 0.0017 


1.42 + 0.14 


0.19 


+ 0.05 


17.2 


F3-Fori, 


F5 


2.3520 ± 0.0022 


0.72 + 0.09 


0.01 


+ 0.06 


13.6 


F1+F2 


F6 


0.0252 ± 0.0032 


0.67 + 0.13 


0.41 


+ 0.09 


9.2 


If 


F7 


2.3964 ± 0.0026 


0.51 + 0.07 


0.48 


+ 0.07 


11.6 


F1+F3 


F8 


2.4486 ± 0.0027 


0.42 + 0.07 


0.24 


+ 0.07 


11.0 


2F1 


F9 


2.3015 ± 0.0026 


0.42 + 0.06 


0.53 


+ 0.07 


11.2 


F2+F3 


FIO 


0.0566 ± 0.0055 


0.45 + 0.14 


0.17 


+ 0.15 


5.4 


F1-F3 


Fll 


1.2117 + 0.0045 


0.36 ± 0.09 


0.67 + 0.12 


6.6 


Fl 


F12 


0.0409 + 0.0050 


0.28 ± 0.08 


0.97 + 0.13 


5.9 


If 


F13 


2.1709 + 0.0035 


0.30 + 0.06 


0.49 


+ 0.09 


8.5 


Fl+F3-F„,b 


F14 


0.4666 + 0.0063 


0.29 + 0.10 


0.96 


+ 0.17 


4.8 


2 Forb 


F15 


1.1133 + 0.0047 


0.26 ± 0.07 


0.65 


+ 0.13 


6.4 


2F3-F1 


F16 


2.0726 + 0.0043 


0.25 + 0.06 


0.77 


+ 0.11 


7.0 


F2+F3-F„,b 


F17 


2.1161 +0.0041 


0.23 + 0.05 


0.52 


+ 0.11 


7.3 


Fl+F2-Fo,b 


F18 


1.1925 + 0.0055 


0.21 + 0.07 


0.36 


+ 0.15 


5.4 


F3-F6 


F20 


1.2744 + 0.0061 


0.16 + 0.06 


0.69 


+ 0.16 


4.9 


2F1-F3 


F21 


2.2493 + 0.0045 


0.15 + 0.04 


0.67 


+ 0.12 


6.6 


2F2 


F22 


3.5254 + 0.0041 


0.15 + 0.03 


0.27 


+ 0.11 


7.3 


F1+F2+F3 


F23 


0.0113 + 0.0046 


0.24 + 0.06 


0.20 


+ 0.12 


6.5 


If 


F24 


3.5750 + 0.0043 


0.13 + 0.03 


0.03 


+ 0.11 


7.0 


F1+2F3 


F25 


3.4758 + 0.0044 


0.12 + 0.03 


0.54 


+ 0.12 


6.7 


F2+2F3 


F26 


1.1394 + 0.0065 


0.15 + 0.06 


0.72 


+ 0.18 


4.6 


F1-2F2-2F3 


F27 


3.6220 + 0.0047 


0.10 + 0.03 


0.64 


+ 0.13 


6.3 


2F1+F3 


F28 


3.2991 + 0.0060 


0.09 + 0.03 


0.60 


+ 0.16 


5.0 


Fl+F2+F3-Fo,b 



Closest frequency combination. The independent frequencies are denoted with "g", the low frequency component with "If. 



stellar rotation and the beating is due to slightly different rotation 
rates of migrating spots. 

In our case, however, even if the oscillatory pattern does 
show symmetry with respect to its mean value, we can exclude 
this hypothesis, because spectroscopy allowed an accurate mea- 
surement of Vrot sin i (and we have from the light curve analysis 
the value of the inclination). Both stars have rotation frequen- 
cies close to the orbital one (1.1 and 0.9 Forb for the primary and 
secondary, respectively), and Forb - Frot is much smaller than 
the dominant frequency of the amplitude spectrum. If the rota- 
tion axes were not aligned with the orbital one, the stellar rota- 
tion frequencies could be higher, but the disalignement should 
be conspicuous: the factor ^ 5 between the dominant frequen- 
cies and Forb requires an inclination of the rotation axis as low 
as 10°. Besides, the similarity between the rotation frequencies 
- as derived in the hypothesis of alignment - and the orbital fre- 
quency would remain an unexplained coincidence. 

We adopt, therefore, the simplest scenario in which the pe- 
riodic variations are not due to rotation but to stellar pulsations. 
The dominant frequencies seen in the spectrum (Fl - F3) fall 
indeed in the range of excited high order g-modes typical of y 
Dor pulsators. 

W e have, therefore, computed with MAD (^ Dupret et al.l 
l2003h the instability domain for stellar models with the physi- 
cal characteristics of the components. In particular, in the case 
of models without overshooting the range is 0.50 - 4.16 (/' = 1) 
and 0.52 - 5.52 d^H^ = 2) for the piimaiy, and 0.55-4.7 (t = 1) 
and 0.55 - 5.42 d"'(^ = 2) for the secondary component. The 
values change only slightly for the models with overshooting, 
becoming 0.51 - 3.83 {I = 1) and 0.55 - 5.02 dr^(t = 2) for the 
primary, and 0.55 - 4.16 (i = 1) and 0.53 - 5.33 d'\{ = 2) for 
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Fig. 10. The standard deviation, with respect to the local mean 
of the binned residual light curve (see text). The vertical lines 
bound the in-eclipse phases. The dip around phase zero suggests 
that the pulsations belong to the eclipsed component, i.e. to the 
primary star 



the companion. Being the range very similar for both stars it is 
impossible to identify the pulsator in this way. 

We have found, however, a possible indirect evidence that the 
pulsating star is the primary, or at least that the higher amplitude 
frequencies belong to the primary. This might be, in fact, the 
reason why the dispersion of the phased light curve decreases 
during eclipse. To show this effect we re-binned the residual light 
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curve (after subtraction of the best binary model) in 200 equal 
phase-bins and computed the local standard deviation, which is 
shown in Fig. [10] as function of the orbital phase. There is an 
evident dip around phase zero. However, taking into account that 
we have only thirteen (narrow) minima in the whole light curve, 
implying a poor sampling of pulsation in eclipse, and that only 
a small fraction of the primary surface is out of view during a 
grazing eclipse, we think that this conclusion shall be taken with 
caution. 

Most frequencies shown in Table|6]can be re-conducted to a 
combination of the first three ones and the orbital frequency. The 
main features in the spectrum are three groups centered around 
F3, F5 and F22 (1.1716, 2.3520 and 3.5254 d"', respectively). 
Each of those frequencies is the center of a structure with side 
peaks spaced by about 0.05 d"' . This fact could suggest an inter- 
pretation in terms of rotational splitting of the central frequency. 
In particular both F5 and F22 have two peaks each side (one 
in the case of F22 just below the significance threshold, Fig|9]). 
This occurence could be interpreted - at first sight - as quin- 
tuplets due the effect of rotation of w\ t - 2 mode, which are 
considered the ones preferen tially excited by tidal forces (e.g., 
iKosovichev & Severnviill983h . 

There are, however, several arguments against this hypothe- 
sis. The first is, again, the known values of the component rota- 
tion. As the star is rotating slowly enough (its rotation frequency 
is 1.5 % of the critical value as defined by Townsend (2005) and 
smaller than the dominant pulsation frequencies), the splitting 
can be express ed according to the first order approximation by 
lLedouxl(fl95l : 



where I is the degree of the mode and F^^i - 0.253 d 'the rota- 
tion frequency. In our case, with 1-2, KF ^ 0.21 d"', i.e. about 
four times the measured value. In the case of a dipolar mode, 
I - \, the resulting splitting of the corresponding triplet would 
be smaller but still ~2.5 times larger than the derived value. 

Even if we interpret the measured value as reflecting the av- 
erage value of the stellar rotation (which is probably not rigid), 
we would end up with a star whose surface (observed) rota- 
tion rate is much faster than that of its interior That might 
be the case of late-type (solar and later) components in close 
binaries, which might be spun up starting from a slow rota- 
tion regime as a consequence of angular momentum loss by 
magnetic braking and spin-orbit synchronization by tidal ac- 
tion. Depending on the core-surface coupling, the external layers 
could, in this case, rotate faster than the interior The components 
of CoRoT 102918586 , however, are early F stars, so the previ- 
ous scenario does not apply. 

Besides, it has to be noticed that the actual value of the spac- 
ings F3-F2 and F1-F3 are not exactly 0.05 but rather 0.0461 + 
0.0016 and 0.0531 ± 0.0020, respectively and, again on the basis 
of the slow rotation, one would expect symmetrical peaks. 

Last, but not least, there exists an alternative and more con- 
vincing interpretation of the central frequencies as combina- 
tions: F5 is closer to F1+F2 than to 2F3 and F22 to F1+F2+F3. 
The remarks in Table |6] correspond to the best matching combi- 
nations. So we can safely conclude that the patterns in Fig.|9]are 
not related to stellar rotation. 

The comparison of the frequency spectrum with those of 
other y Dor stars, observed by CoRoT and Kepler, suggests that 
the pattern we see is quite common in (pre sumably ) single stars 
as well (see, e.g., iBalona et al]l201 iL iHareter et al.ll2010l) . For 



instance, the light curv^ of CoRoT 102732872 (an F3V star, 
pulsating with a main frequency of 1.09 d"') closely resembles 
that of CoRoT 102918586 after removal of the EE contribution. 
The former star is one of the few j Dor which, according to 
IHareter et al. shows a pattern of regularly spaced peaks 

in period (as expected in the case of high-order g-modes). We 
computed, therefore, the mean period spacing of the three inde- 
pendent frequencies, AP = O!0360 ± OJOOl 1 (or 3 1 10'^ + 900 and 
compared it with the theoretical asymptotic period spacing of the 
i - \ high-order g-modes, as expected for the stellar parameters 
from Table[3]and Section|6] 

The computed range of period spacing for the primary star, 
relative to the stellar models in the age interval derived in 
Section |6] is 3070-3200^ for the models without overshooting, 
and 3300-3390^ for those with aov - 0.2. For the secondary 
star the values are, respectively 2800- 3030' and 3160 - 3300'. 
A consistent interpretation is, therefore, that the observed pattern 
is due to non-radial dipolar pulsations of the primary component. 

8. Discussion and conclusions 

CoRoT 102918586 is the first y Dor in an eclipsing binary for 
which it was possible to collect observations of such a quality 
to allow a very detailed study, which provided both an accurate 
determination of the physical parameters and of the pulsational 
properties. Masses and radii have uncertainties < 1%, the effec- 
tive temperatures of 1.5 % and the derived age, turned out to be 
600 Myr, within a range of ~ 100 Myr 

Our results depict the system as a young eccentric binary 
formed by two similar MS star (the mass ratio is 0.9), the sin- 
gle eclipse is just due to the orientation of the binary orbit in 
space and the primary component is at the origin of the observed 
y Dor pulsations. We cannot completely exclude that the similar 
secondary component pulsates as well, but we suggest the dom- 
inant pulsation frequencies belong to the primary, on the basis 
of the measured decrease of the pulsation amplitude during the 
eclipse. 

Spectra disentangling of the FEROS spectra allowed to 
derive the effective temperatures and the component chemi- 
cal abundances, which turned out to be close to solar The 
abundance pattern for the y Dor pulsating stars was studied 
by Br untt et al. (2008), and more recently by Tkache nko et al.l 
( 120121) . particularly in regard to the links with t he chemically 
peculi ar A Boo and Am-type stars suggested by iGrav & Kavd 
(Il999h and ISadakanj (l2006l) . Their comprehensive detailed 
abundance analysis did not reveal any sign of chemical peculiar- 
ities in their sample of y Dor stars, yielding the conclusion that 
the metallicity is quite close to the solar value. P reviously , on the 
basis of the Stromgren photometric m\ index, iHandled (Il999h 
as well noted that all the y Dor candidates in his sample had 
metallicities close to the solar value. Our finding for the y Dor 
pulsating component in the binary system CoRoT 102918586 
corroborates this general conclusion. 

The results of our asteroseismic study are limited by the lack 
of calibrated multicolor photometry in CoRoT and/or of spec- 
troscopy of sufficiently high S/N to allow to study line profile 
variations due to pulsation (unfeasible given the system magni- 
tude). As a consequence, we could not completely identify the 
pulsation modes. On the other hand, from the comparison of the 
period splitting with the theoretically expected values, we have 

- The light curve and its plot are available from the public sec- 
tion of the CoRoT-N2 archive at the CoRoT Data Center (IAS): 
http://idoc-corot.ias.u-psud.fr/, 
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at least indications that the primary component is pulsating in a 
I - \ mode. Those are not the modes preferentially excited by 
tidal f orces, which h a ve bee n detected (or suggested) in in other 
cases. [Handler et al ] (120021) propose that a number of frequen- 
cies, exact multiple of the orbital frequency, found in the y Dor 
component of HD 209295 (a single-lined spectroscopic binary) 
could be triggered by tidal interaction. Another case is that of 
HD 174884 (Maceroni et al. 2009), an eccentric eclipsing binary 
formed by two B-type stars. 

In CoRoT 102918586 , however, we find little or no influence 
of binarity on the pulsations: the pulsation pattern is very similar 
to that of single y Dor's and the only sign of the orbital motion is 
found in the beating of the orbital frequency with the dominant 
ones. This might be due to a weak tidal interaction: our system 
is less eccentric than the above-mentioned ones, has a longer 
period, and the components have smaller fractional radii (by a 
factor of ~ 1/2 with respect to HD 174884 and even smaller in 
the case of HD 209295, whose components radii - however - 
can only be estimated). One expects, therefore, a weaker tidal 
torque, as this is very strongly dependent on the star fractional 
radius (e.g., Zahn 2005). 

The final model of the system, characterized by an internal 
consistency among the results from diff'erent data set and treat- 
ments, proves the success of the iterative method used to sepa- 
rate pulsations from the eclipsing binary light curve. That was 
not obvious, and we do not think that this result is generally 
valid, because eclipses certainly modify the observed oscillatory 
pattern, while we subtracted the same oscillations at all phases. 
A reason of the success in this case is that we dealt with a graz- 
ing eclipse, which decreases the determinacy of the light curve 
solution but makes this particular problem less severe. 

A straightforward, first-order, correction could be to weight 
the pulsation amplitude with the contribution to the total light of 
the pulsating component during eclipses, but the definite solu- 
tion is the inclusion of non-radial pulsation in the binary model, 
taking into account the deformation of the surfaces by non-radial 
pulsations and the consequent brightness variation. That is how- 
ever a challenging task, as we do not know a priori in which 
modes the star oscillates, and adding the parameters governing 
the oscillations to the minimization procedure greatly increases 
its complexity (and its unicity problems). 
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